Measurement of the double-/? decay half-life of 136 Xe with the KamLAND-Zen experiment 
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We present results from the KamLAND-Zen double-beta decay experiment based on an exposure of 77.6 days 
with 129 kg of 136 Xe. The measured two-neutrino double-beta decay half-life of 136 Xe is Ty 2 = 2.38 ± 
0.02(stat) ± 0.14(syst) x 10 1 yr, consistent with a recent measurement by EXO-200. We also obtain a lower 
limit for the neutrinoless double-beta decay half-life, T°Y 2 > 5.7 x 10 24 yr at 90% confidence level (C.L.), 
which corresponds to almost a five-fold improvement over previous limits. 



PACS numbers: 23.40.-s, 21.10.Tg, 14.60.Pq, 27.60.+j 
I. INTRODUCTION 

Majorana neutrinos are a natural feature of many high- 
energy physics theoretical models. However, the only viable 
experimental probe of this property at present is neutrinoless 
double-beta (0v(3(3) decay Jit]. Observation of this lepton- 
number violating nuclear process would definitively establish 
the Majorana nature of the neutrino, and would be a pro- 
found discovery 0]. In addition, since the rate of this pro- 
cess increases with the square of the effective neutrino mass 
{ m pp) = |^i^ei m ^i|' lts measurement would provide in- 
formation on the absolute neutrino mass scale. Searches for 
Ov(3(3 decay have been invigorated by recent measurements 
of neutrino mass splittings by oscillation experiments, which 
require at least one neutrino mass above ^50 meV Hfl. This 
scale is within the reach of present-day efforts. 

Determining (mpp) from a Qv(3(3 decay half-life requires 
knowledge of the decay's phase-space factor (G 0l/ ) and nu- 
clear matrix element (M°"). G 0l/ can be calculated exactly, 
but to date all estimations of M 0v must rely on model-based 
approximations with difficult-to-quantify uncertainties. The 
two-neutrino double-beta {2v(3(3) decay half-life, if known, 
can be used to constrain some relevant model parameters, 
reducing some sources of uncertainty 0,111]. The first di- 
rect measurement of the 2v(3J3 decay half-life of 136 Xe, re- 
cently reported by EXO-200 |aL_ was significantly below pre- 
viously published lower limits |0, H] ■ This article on the first 
results from the KamLAND-Zen (KamLAND Zero-Neutrino 
Double-Beta Decay) experiment reports a new measurement 



of the 136 Xe 2v(3(3 decay half-life, as well as improved limits 
on the 0i/(3(3 mode. The data presented were collected be- 
tween October 12, 201 1, and January 2, 2012. 



II. DETECTOR AND CALIBRATION 

KamLAND-Zen (Fig. [TJ is a modification of the exist- 
ing KamLAND detector carried out in the summer of 2011. 
The (3(3 source/detector is 13 tons of Xe-loaded liquid scin- 
tillator (Xe-LS) contained in a 3.08-m-diameter spherical in- 
ner balloon (IB). The IB is constructed from 25-/im-thick 
transparent nylon film and is suspended at the center of the 
KamLAND detector U by 12 film straps of the same material. 
The IB is surrounded by 1 kton of liquid scintillator (LS) con- 
tained in a 13-m-diameter spherical outer balloon (OB) made 
of 135-/xm-thicknylon/EVOH (ethylene vinyl alcohol copoly- 
mer) composite film. The outer LS is 0.10% less dense than 
the Xe-LS and acts as an active shield for external 7's and as a 
detector for internal radiation from the Xe-LS or IB. The Xe- 
LS consists of 82% decane and 18% pseudocumene (1,2,4- 
trimethylbenzene) by volume, 2.7 g/liter of the fluor PPO (2,5- 
diphenyloxazole), and (2.52 ± 0.07) % by weight of enriched 
xenon gas, as measured by gas chromatography. The isotopic 
abundances in the enriched xenon were measured by residual 
gas analyzer to be (90.93 ±0.05)% 136 Xe and (8.89 ±0.01)% 
134 Xe; other xenon isotopes are negligible. The light yield of 
the Xe-LS is 3% lower than that of the LS. Buffer oil (BO) be- 
tween the OB and an 18-m-diameter spherical stainless-steel 
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FIG. 1 : Schematic diagram of the KamLAND-Zen detector. 



containment tank (SST) shields the LS from external radia- 
tion. Scintillation light is recorded by 1,325 17-inch and 554 
20-inch photomultiplier tubes (PMTs) mounted on the SST, 
providing 34% solid-angle coverage. The SST is surrounded 
by a 3.2 kton water-Cherenkov outer detector (OD). Details of 
the KamLAND detector are given in Ref . ^ . 

The data acquisition system (DAQ) is triggered when 70 
or more 17-inch PMTs are hit (primary trigger), which cor- 
responds to a threshold of ~0.4 MeV. The signals on all hit 
PMTs are digitized for ^200 ns for offline analysis. After 
each primary trigger the threshold is lowered to ~0.25 MeV 
for 1 ms to study sequential decays. The scintillation light 
from the two coincident e~ produced by 136 Xe (3f3 decay can- 
not be separated, so only their summed energy is observed. 
For hypothetical 0^/3/5 decays, the sum is always 2.458 MeV 
(Q value of the 136 Xe /3/3 decay) d, while for the 2^/3 
decays the sum has a continuous spectrum up to the Q value. 
Event energy (visible energy) is estimated from the number 
of observed photoelectrons (p.e.) after correcting for PMT 
gain variation and solid angle, shadowing, and transparency 
of detector materials. The corrections depend on the event 
vertex. The vertex reconstruction is based on the maximum 
likelihood fit to the pulse shape of each PMT hit timing after 
correcting for photon time of flight. The pulse shape is al- 
most determined by scintillation decay time and dark hit con- 
tribution, and it differs between 17-inch and 20-inch PMTs 
due to different transit-time spreads. The vertex resolution 
is estimated from radial distributions of radioactiv e contam- 
inants (see Fig. to be a ^15 cm/ \J ' E{McV). The en- 
ergy response is calibrated with 7's from a 208 T1 (TI1O2W) 
source, 214 Bi (/? + 7's) from 222 Rn (t = 5.5 day) introduced 
during detector modification, and 2.225 MeV 7's from spalla- 
tion neutrons captured on protons. 

Figure [2 a) shows the energy spectrum obtained when the 
TI1O2W source, contained in a ^5-mm-thick lead capsule, 
was deployed close to the outer surface of the IB. The most 
intense peak is due to the primary 7 of 208 T1 (2.614 MeV). 
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FIG. 2: Visible energy distributions of (a) 7's from the 208 T1 cal- 
ibration source and (b) 214 Bi (/3 + 7's) decays in the Xe-LS. The 
lines indicate the best fits to the analytical spectral models with the 
resolution and energy scale parameters floating. The fit to (a) has a 
X 2 /d.o.f. = 5.0/8 and (b) has a x 2 /d.o.f. = 27.0/29. 



The less intense peak near ^3.1 MeV is from multiple-7 cas- 
cades of 208 T1. According to Monte Carlo (MC) studies, the 
degradation of the primary 7 inside the source is negligible, 
and the distribution around the primary peak can be described 
by a Gaussian distribution and a third-order polynomial. The 
mean and width of the Gaussian distribution are relatively in- 
sensitive to the polynomial parameters. The resultant energy 
resolution at 2.614 MeV is a = (6.6 ± 0.3)%/ y/E(MeV). 
The parameters of a detector energy nonlinear response model 
describing effects from scintillator quenching and Cherenkov 
light production are constrained to reproduce the 2.614 MeV 
208 T1 peak position and the spectral shape of 214 Bi events 
[Fig- 0t))]. From the neutron-capture 7 data, systematic vari- 
ation of the energy reconstruction over the Xe-LS volume is 
less than 1.0%, and the detector energy response is stable to 
within 1.0% during the data set. 



III. CANDIDATE EVENT SELECTION 

Candidate /3j3 decay events are selected by performing the 
following series of cuts: (i) The reconstructed vertex must be 
within 1 .2 m of the detector center, defining the fiducial vol- 
ume (FV). (ii) Muon (events with more than 10,000 p.e. or 
more than 5 OD hits) and events occurring within 2 ms af- 
ter muons are eliminated, (iii) A coincidence cut eliminates 
sequential events that occur within 3 ms of each other; this re- 
moves (99.97 ± 0.01)% of 214 Bi- 214 Po (/3 + 7, then a decay, 
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r=237 fis) decays. This cut is augmented with a secondary 
cut, aimed at identifying sequential 212 Bi- 212 Po (t=0.4 /is), 
which exploits detailed PMT wave-form data to identify co- 
incidences within a single ^200-ns-long DAQ event window. 
The 212 Bi- 212 Po rejection efficiency is (89 ± 2)%. The vis- 
ible energy of a decay is quenched due to the high ioniza- 
tion density, so the cuts are applied only if the delayed a's 
are between 0.35 and 1.5 MeV. The dead time introduced by 
the coincidence cuts is less than ~0.1%. (iv) A background 
mainly from reactor 77 e 's producing a delayed coincidence of 
positrons and neutron-capture 7's (2.225 MeV) is rejected by 
requiring event time separations greater than 1 ms, and the de- 
layed energy is larger than 1.5 MeV. (v) Finally, candidates 
must pass a vertex-time-charge (VTQ) test designed to fil- 
ter out noise events. The test compares the observed PMT 
charge and hit-time distributions to those expected based on 
the reconstructed vertex full . The VTQ cut is tuned using 
KamLAND calibration data and reduces the selection effi- 
ciency by less than 0.1%. The total livetime after all cuts is 
77.6 days. The energy spectrum of (3(3 decay candidates is 
shown in Fig. [4] 



IV. BACKGROUND ESTIMATION 

Backgrounds to the /3/3 decay study fall into three cat- 
egories: those external to the Xe-LS, mainly from the IB 
material; those from residual radioactive impurities in the 
Xe-LS; and spallation backgrounds (Table [Q. From a spec- 
tral analysis of events whose reconstruction is close to the 
IB boundary, we find that the activity in the energy region 
1.2 MeV < E < 2.0 MeV (2^/3/3 window) is dominated by 
134 Cs (0 + 7); in the region 2.2 MeV < E < 3.0 MeV (0i//3/? 
window), the spectrum is consistent with 214 Bi (/3+7) decays. 
The observed surface activity ratio of 134 Cs to 137 Cs (0.662 
MeV 7) is ~0.8 which is consistent with contamination by 
fallout from the Fukushima-I reactor accident in March 201 1. 
The IB production facility is located just 100 km away from 
the Fukushima-I reactor. The FV cut is performed to miti- 
gate the background from the IB material; the remaining IB 
background inside the FV is estimated by fitting Monte-Carlo- 
generated event radial distributions to the data. Figure[3]shows 
the event density as a function of the cubed radius from the IB 
center for the two energy ranges, along with fits to the MC 
distributions. In the 2z//3/3 window we fit for a 2vf3j3 source 
uniformly distributed in the Xe-LS and a 134 Cs source uni- 
formly distributed on the IB. In the 0^/3/3 window we show 
the best fits for a 214 Bi source uniformly distributed on the 
IB and either a Of /3/3-like source or a 2.6 MeV 7 source uni- 
formly distributed in the Xe-LS. The radial distribution offers 
no discrimination between these event types. 

Assuming secular equilibrium, the residual 238 U and 232 Th 
concentrations internal to the Xe-LS are estimated to be 
(3.5 ± 0.6) x lO" 16 g/g and (2.2 ± 0.3) x 10~ 15 g/g, respec- 
tively, based on sequential decays of 214 Bi- 214 Po and 212 Bi- 
212 Po. Since equilibrium may be broken by introduction of 
contaminants during detector modification, the Bi-Po stud- 
ies are only used to estimate internal background from the 



TABLE I: Summary of background in the fiducial volume. Note that 
134 Xe /3/3 contribution is expected to be negligible because of much 
smaller phase space factor and isotopic abundance relative to 136 Xe. 
The backgrounds from nuclear reactions such as (a, 7) and (a, cvy) 
are expected to be negligible due to their small cross sections. The 
rate of (n, 7) is stringently constrained from the data. 



Isotope 



Event rate (ton day) 



External (radioactivity in IB) 



238-j 

232 Th series 22S Th- 208 Pb 

40 K 

4 Cs 
r Cs 



134( 
137,- 



(3.2 ±0.3) x 1(T 2 
(3.5 ±0.3) x 10~ 2 
0.21 ±0.16 
0.50 ±0.07 
0.35 ±0.06 



Residual radioactivity in Xe-LS 



5T T 

U series 



238 TJ- 222 Rn ( 234 Pa) 



232- 



2 Rn- 210 Pb 

228 



Th series 



Th- 22S Th ( 22S Ac) 



^Th- 208 Pb 



'Kr 



228^ 208 T 

85 j 
210 Bi 
40 K 

Cs 
Y 
Cs 



134 
90 
137 



<1.5 
4.9 ± 0.2 

<0.7 
0.58 ± 0.06 
196 ±8 
103 ±3 

<9.6 

<0.4 

<0.8 

<1.1 



Spallation product from 12 C a 



10 C (2.11 ±0.44) x 10~ 2 

U C 1.11 ±0.28 

Spallation products from xenon with lifetime < 100 s 



1.2 MeV < E < 2.0 MeV 
2.2 MeV < E < 3.0 MeV 



<0.3 
<0.02 



Around 136 Xe Ou/3/3 (ENSDF search) 



60 Co, 88 Y, 110 Ag m and 208 Bi 



0.22 ±0.04 

Potential background from fallout 

listed above 
negligible 



134 Cs, 137 Cs, 110 Ag m 

129 j g m 95-ntk 90 v 89 1 



, Nb «U Y «« Sr 



a Other spallation products are negligible jjj. 



222 Rn .2io pb su bchain of the 238 U series, and from the 228 Th- 
208 Pb subchain of the 232 Th series; other decays in both series 
are treated as unconstrained backgrounds. We note the well- 
known 2.614 MeV 7 from 208 T1 (/3~ decay, Q = 5.00 MeV) 
in the 232 Th series is not a serious background for the Oi//3/3 
decay search because of detection of the coincident (3/j in the 
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FIG. 3: R 3 vertex distribution of candidate events for (a) 1.2 MeV < 
E < 2.0 MeV and (b) 2.2 MeV < E < 3.0 MeV. The curves 
show the best-fit model components: (a) 2v/3/3 (dashed) and 134 Cs 
(dotted), (b) 2.6 MeV 7's (dashed) and 214 Bi (dotted); the long- 
dashed curve is for Ou/3/3 instead of 7's. (c) 214 Bi events from Xe-LS 
(dashed) and IB (dotted). The vertical lines show the fiducial radius 
of 1.2 m (dashed) and the IB radius (dotted). 



surrounding active LS 1L 

Spallation neutrons are tagged by coincidence of neutron- 
capture 7's with preceding muons. We expect capture on pro- 
tons (2.225 MeV), 12 C (4.946 MeV), 136 Xe (4.026 MeV), 
and 134 Xe (6.364 MeV), with fractions 0.994, 0.006, 9.5 x 
10~ 4 , and 9.4 x 10~ 5 , respectively. We find no 136 Xe or 
134 Xe neutron-capture candidates in the data set. 137 Xe (/3~, 
t == 5.5 min, Q = 4.17 MeV) from neutron capture on 136 Xe 
is a potential 0i//3j3 background, but the expected production 
rate is negligible, ~2.9 x 10~ 3 (ton day) -1 , where ton is a 
unit of Xe-LS mass. Production rates of light nuclei by spal- 
lation of carbon are calculated from spallation yields previ- 
ously measured in KamLAND U. We observe a (13 ± 6)% 
increase in the spallation neutron flux in the Xe-LS rela- 
tive to the outer LS, from which we assess a 19% system- 
atic uncertainty on the calculated spallation yields. Events 
from decays of n C (/3+, r = 29 A min, Q = 1.98 MeV) and 
10 C (/3+, t = 27.8 s, Q = 3.65 MeV) dominate the contri- 
butions from spallation backgrounds. We expect rates of 
l.lliO^Sttonday)- 1 and (2.11±0.44) x 10" 2 (ton day)" 1 
from 11 C and 10 C, respectively. The 11 C/ 10 C background 
can be reduced by a triple-coincidence tag of a muon, a neu- 
tron, and the 11 C/ 10 C decay. This is not pursued in the cur- 
rent analysis. We found no past experimental data for muon 
spallation of xenon. With the present data, we find the event 
rates in the energy ranges 1.2 MeV < E < 2.0 MeV and 
2.2 MeV < E < 3.0 MeV from isotopes with lifetimes of 
less than 100 s associated with muons depositing more than 
~3 GeV in the detector (so-called showering muons) are less 
than 0.3 (ton day) _1 and 0.02 (ton day)" 1 at 90% C.L., re- 



spectively. 

Care has to be taken for backgrounds that produce a peak 
close to the Oz//3/3 decay energy (2.458 MeV), particularly 
ones which may have been introduced during detector mod- 
ification or may be induced by muon spallation. We searched 
all isotopes in the ENSDF database lfl3ll for sources with 
a peak structure between 2.4 and 2.8 visible-MeV. For all 
tabulated decay chains we calculate the visible energy spec- 
trum, accounting for the time structure of the chain and pileup 
in the DAQ event window, and apply the nonlinear detec- 
tor energy response model to all individual decay secon- 
daries of each branch. Considering only nuclei with decay 
ancestor lifetimes longer than 30 days, we identify 110 Ag' n 
(/?- decay, r = 360 day, Q = 3.01 MeV), 88 Y (electron cap- 
ture (EC) decay, r = 154 day, Q = 3.62 MeV), 208 Bi (EC de- 
cay, r = 5.31 x 10 5 yr, Q = 2.88 MeV), and 60 Co 0~ de- 
cay, r = 7.61 yr, Q = 2.82 MeV) as potential background 
sources. Observation of 134 Cs/ 137 Cs on the IB raises the plau- 
sibility of contamination of detector materials by Fukushima 
fallout, which includes 110 Ag m . One assay of soil samples 
taken near the IB production facility revealed evidence of 
no Ag m Although 88 Y, 208 Bi, and 60 Co are not detected near 
Fukushima or in our soil samples, we consider them to be pos- 
sible backgrounds. Except for 208 Bi, these long-lived back- 
ground candidates can be also produced from xenon spalla- 
tion by cosmic -rays when materials were aboveground, but 
the rate estimations are difficult. Broadening the search to in- 
clude shorter-lived nuclei (100 s < r < 30 day) possibly 
supported by muon spallation in the detector, we found that 
the production of candidate parents with mass numbers be- 
low 136 Xe is stringently constrained by comparing production 
cross sections in Ref. IU4I1 . 



V. SYSTEMATIC UNCERTAINTY 



Nominally, the 1.2-m-radius FV corresponds to 0.438 ± 
0.005 of the total Xe-LS volume (16.51 ± 0.17 m 3 ), or 129 kg 
of 136 Xe. The fiducial volume fraction may also be esti- 
mated from the fraction of 2l4 Bi events whose reconstruction 
is within 1 .2 m of the IB center compared to the total number 
in the entire Xe-LS volume after subtraction of the IB surface 
contribution. The result is 0.423 ± 0.007(stat) ± 0.004(syst) 
[Fig-EJc)], including the uncertainty on the IB surface contri- 
bution. The slight event loss around the center may indicate an 
outward radial vertex bias, so the difference in these estimates 
is taken as a measure of the systematic error on the vertex- 
defined FV. Combining the errors, we obtain a 5.2% system- 
atic error on the fiducial volume. The total systematic uncer- 
tainty on the /3/3 decay half-life measurement is 5.9%, com- 
ing from the quadrature sum of the fiducial volume (5.2%), 
enrichment of 136 Xe (0.05%), Xe concentration (2.8%), de- 
tector energy scale (0.3%), Xe-LS edge effect (0.06%), and 
detection efficiency (0.2%). 
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FIG. 4: (Color) (a) Energy spectrum of selected /3/3 decay can- 
didates together with the best-fit backgrounds and 2v/3/3 decays, 
and the 90% C.L. upper limit for 0^/3/3 decays; the fit range is 
0.5 MeV < E < 4.8 MeV. The coincident /3/~/ events from 
208 T1 give the main contribution from 3 to 5 MeV. (b) Closeup 
of (a) for 2.2 MeV < E < 3.0 MeV. (c)-(e) Closeups for 
2.2 MeV < E < 3.0 MeV with different background assumptions 
around the OvfiP peak: (c) only OufiP, (d) 110 Ag m + 0^/3/3, (e) 
2c,8 Bi+ Of /3/3, for comparison, 
range. 
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VI. RESULT 

The 136 Xe 2isf3(3 and Of /3/3 decay rates are estimated from 
a likelihood fit to the binned energy spectrum of the selected 
candidates between 0.5 and 4.8 MeV. The 136 Xe 2v/3P de- 
cay spectrum shape from Ref. lfl5tl is used; 80% of the 2f /3/3 
spectrum falls within our energy window. The contributions 
from major backgrounds in the Xe-LS, such as 85 Kr, 40 K, 
nonequilibrium 210 Bi, and the 238 U- 222 Rn and 232 Th- 228 Th 
decay chains, are free parameters and are left unconstrained 
in the fit. The contributions from the 222 Rn- 210 Pb and 228 Th- 
208 Pb chains, 11 C, and 10 C are allowed to vary but are con- 
strained by their independent measurements. Residual IB- 
surface backgrounds in the FV are constrained by the ra- 
dial distribution study. The parameters of the detector en- 
ergy response model are floated but are constrained to re- 
produce the 208 T1 source and 222 Rn-induced 214 Bi data. Po- 
tential backgrounds from fallout nuclei with half-lives longer 
than 30 days found in ex situ measurements of soil or ocean 
samples around Fukushima, namely 137 Cs, 134 Cs, 110 Ag m , 
129 Te m , 95 Nb, 90 Y (from 90 Sr), and 89 Sr, as well as potential 
Oz//3/3 backgrounds found in the ENSDF search ( 88 Y, 208 Bi, 
and 60 Co) are included as unconstrained free parameters. The 
relative contributions of Of /3/3 window backgrounds are addi- 
tionally constrained by the time variation of the event rate in 
the energy range 2.2 MeV < E < 3.0 MeV. 

Figured! a) shows the resulting best-fit spectral decomposi- 
tion. The x 2 /d.o.f. comparing the binned data and the best-fit 
expectation is 99.7/87. 2f/3/3 decay is the dominant spec- 
tral feature in the low-energy region. The best-fit number 
of 136 Xe 2f/3/3 decays is (3.55 ± 0.03) x 10 4 , correspond- 
ing to an event rate of 80.9 ± 0.7 (ton day) -1 . We found 
no systematic variations due to the choice of the data period 
and volume within the 1 .2-m-radius FV. The dominant back- 
grounds at low-energy are from 85 Kr and 210 Bi, with best- 
fit rates of 196 ± 8 (ton day) _1 and 103 ± 3 (ton day)" 1 , re- 
spectively. The fit yields the following 90% C.L. upper limits 
on other background rates (per ton-day) in the Xe-LS: 40 K 
<9.6, 234 Pa<1.5, 134 Cs<0.4, 228 Ac <0.7, 90 Y < 0.8, and 
137 Cs<l.l; other fallout isotopes are negligible. 

Around the Of/3/3 energy, a strong peak appears, but the 
peak is centered significantly above the Q value of the de- 
cay [Fig.Uc)]; the hypothesis that the peak can be described 
by Of/3/3 decay alone is rejected by a \ 2 test at more than 
5(7 C.L., including the systematic uncertainties on the energy 
scale model. The best-fit combined background rate around 
the 0^/3/3 energy allowing for contributions from 110 Ag m , 
88 Y 208 Bi, and 60 Co is 0.22 ± 0.04 (ton day) -1 . Figures|ld) 
and 2te) show the distribution if only 208 Bi or 110 Ag" 1 , re- 
spectively, contribute as background to the Of /3/3 peak. We 
conclude that the data in the Of/3/3 region is contaminated 
by the candidate backgrounds (but mainly due to 208 Bi or 
110 Ag m ), and the Of/3/3 limit is extracted by floating those 
contributions [Fig. Htb)]. The 90% C.L upper limit on the 
number of 136 Xe Of/3/3 decays is < 15 events, an event rate 
of < 0.034 (ton day)" 1 . 

The measured 2f/3/3 decay half-life of 136 Xe is T\» = 

2.38 ± 0.02(stat) ± 0.14(syst) x 10 21 yr. This is consis- 
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tent with the result obtained by EXO-200, T 2 ^ = 2.11 ± 

0.04(stat) ± 0.21(syst) x 10 21 yr H. For Qvfifi decay, 
the data give a lower limit of T°J 2 > 5.7 x 10 24 yr (90% 
C.L.), which corresponds to almost a fivefold improvement 
over previous limits 01- From the limit on the 136 Xe 0v(3(3 
decay half-life we obtain a 90% C.L. upper limit of (mpp) < 
(0.3 - 0.6) eV using recent QRPA (CCM SRC) d and shell 
model IU7I1 nuclear matrix elements calculated prior to the 
EXO-200 measurement. 



VII. CONCLUSION 

In summary, KamLAND-Zen provides an improved mea- 
surement of the 136 Xe 2i//3/3 decay half-life. The result is 
consistent with that of EXO-200 and supports the conclusion 
that the directly measured half-life is significantly less than 
the lower limits reported by earlier experiments. Our analysis 
includes a search for 0i/(3f3 decay of 136 Xe and yields an im- 



proved lower limit on its half-life. Removal of contaminants 
in the Xe-LS is an important task to improve the 0^/3/3 de- 
cay search sensitivity. In the future, systematic uncertainties 
will also be reduced by performing source calibrations in the 
Xe-LS. 
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